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(57) Abstract: The invention relates to a metiiod 
and device by means of which an irrotational, 
sourceless vector field can be expressed by a 
number of physically reasonable basis vectors 
as small as possible in such a manner that the 
measured signals can be unambiguously divided 
into signals of the irrotational, sourceless vector 
field that are caused by an interesting object or 
external interferences, as well as into a signal 
caused by the nonideality of the measuring device, 
which signal is not included in the model of the 
signal space describing an irrotational, sourceless 
vector field. The invention is based on the 
combining of two very fundamental mathematical 
regularities and applying in the processing of signal 
vectors of a multi-channel measuring device that 
measures an irrotational, sourceless vector field: on 
the Maxell's equations of an irrotational, sourceless 
vector field, as well as on the convei^ence 
characteristics of series developments. 
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Method and system for processing a multi-channnel 
measurement of magnetic fields 

FIELD OF THE INVENTION 

The present invention relates to a new and 
advanced method and system for processing a multi- 
channel measurement signal. In particular, the present 
invention relates to measurement geometry that meas- 
ures a sourceless, irrotational vector field, which 
measurement geometry consists of a sensor arrangement; 
the invention also relates to the formation of a basis 
for signal space and the utilisation of. the basis in 
the pre-processing and analysis of data measured by 
the measurement geometry. 

BACKGROUND OF THE INVENTION 

The representation of measurement signals or 
signal vectors measured by multi- channel measuring de- 
vices in some natural signal space basis is . useful 
from the point of view of the analysis and illustra- 
tion of the measured phenomenon. For example, in the 
case of a multi-channel magnetometer, one such basis 
can be obtained by modeling the source that produced 
the magnetic field by a current distribution whose to- 
tal amplitude is the smallest possible one among the 
distributions that are capable of producing the meas- 
ured signal vector and which current distribution has 
been limited within a conductor area that describes 
the object to be examined. This kind of method is 
called a minimum norm estimate . and it has. been de- 
scribed e.g. in publication "Magnetoencephalography - 
theory, instrumentation, and applications to noninva- 
sive studies of the working human brain" , Matti 
Hamalainen et al Reviews of Modem Physics, Vol. 65, 
No. 2, 1993. In practice, due to the numerical sensi- 
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tivity, the modeling of signal vectors by minimum norm 
estimate requires regularisation whose intensity de- 
pends, for each case specifically, on the signal-to- 
noise ratio and on the distribution of the vector 
field often in a non-obvious manner. Thus, the use of 
minimum norm estimate requires consideration and good 
expertise because a wrong regularisation may lead to a 
completely erroneous modeling. Another problem associ- 
ated with the minimum norm estimate is that it tends 
to model also external interferences in the source 
distribution of an object to be examined, which leads- 
to an erroneous outcome. In addition, the minimum norm 
estimate is a rather laborious operation computation- 
ally. 

The source distribution that produced the 
vector field has also been modelled by means of a 
multi-pole development, but in that case, the source 
model has usually not been formed to be the basis for 
signal space. At a theoretical level, the formation of 
a basis based on a multi-pole development has been 
presented in the doctoral thesis Interpretation of 
Neuromagnetic Measurements: Modeling and Statistical 
Considerations", Matti Hamalainen, 1987, in which the 
basis is designed to describe the magnetic field out- 
side an area containing the source distribution. The 
multi-pole development has also been used for source 
modeling in magnetocardiography, but in that case one 
has not formed a basis for signal space of the devel- 
opment. This is apparent e.g. from the publication 
^'Comparability of Measurement Results Obtained with 
Mult i -SQUID- Systems of Different Sensor Configura- 
tions", M. Burghoff et al, IEEE Transactions on Ap- 
plied Superconductivity, Vol. 7, No. 2, 1997. The for- 
mation of the basis for signal space based on the se- 
ries development of the vector field is reasonable 
only if the measuring device is a sufficiently multi- 
channel one and appropriate with respect to its geome- 
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try to ensure a susceptibility _ to interference as 
small as possible. 

External interferences e.g. from the magnetic 
fields caused by the electric conductors situated near 
the measurement space are easily summed up in subtle 
magnetometer measurements. The elimination of the ex- 
ternal interferences from the measured signals is im- 
portant in order that one could obtain information as 
dependable as possible from the object to be examined. 
Conventionally, as the. interference elimination meth- 
ods, projection and reference signals methods have 
been used. In the former method, the interference 
elimination is based upon information on signal sub- 
space spanned by known interferences, and in the lat- 
ter one upon signals measured by correctly placed ref- 
erence sensors, which signals are assumed to have been 
caused merely by interference sources . 

Of projection methods let it be mentioned the 
SSP projection method in which the sub-space caused by 
typical external interference sources is determined 
and the measured signals are projected against this 
orthogonally to the space. In addition to the inter- 
ference, the projection also diminishes the actual in- 
teresting signal, in case the signal vectors produced 
by the object to be examined are not orthogonal to the 
interference sub-space, and furthermore, the SSP 
eliminates the interference completely only in cases 
in which the external interferences really belong to a 
predetermined interference sub-space. Also the illus- 
tration of the spatial distribution of the signals 
suffers from the projection. The SSP method has been 
described e.g. in patent publication FI925461 and in 
publication '^Signal- space projection method for sepa- 
rating MEG or EEG into components", M.A. Uusitalo and 
R.J. Ilmoniemi, Medical & Biological Engineering & 
Computing, Vol. 35, pp. 135-140, 1997. 
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In reference signal methods, so-called refer- 
ence sensors are installed in the measuring device in 
such places with respect to other sensors that they 
can be considered to measure solely external interfer- 
ences without observing any signal from the object to 
be measured. Taking into account the geometry between 
the reference and the actual signal channels, the sig- 
nals associated with the external interferences can be 
calculated and reduced from the signals measured by 
the signal channels, in which case just the signal re- 
lating to the object to be measured is remaining. This 
kind of method is justified only when the reference 
sensors are capable of measuring all the information 
relating to the external interferences in the area of 
the signal sensors without, however,, measuring the 
signal associated with the object at all. m other 
words, the method is based on an assumption according 
to which exteimal interferences in the area of the 
measuring device are uniform, whereas the signals pro- 
duced by the object to be measured are weakening very 
fast as the distance grows. The method based on the 
reference signals is described e.g. in patent publica- 
tion WO9641209. 

In addition to external interferences, one 
problem that distorts the signals is the possible 
movement of the object to be examined during the meas- 
urement. Lately, e.g. in magnetoencephalography there 
have been made excitation response measurements to 
testees who cannot keep their heads immovable during 
the measurement. An apparent solution is to reject 
from the signals to be averaged signals that corre- 
spond to such moments of time during which the head 
has been situated too far from some reference point. 
In this kind of method, the signal-to-noise ratio 
weakens due to the rejection of the responses particu- 
larly when the head moves much. In more advanced meth- 
ods, the distortion caused by a movement can be cor- 
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rected either by taking the into account the movement 
in the source modeling when using a distorted averaged 
signal, or by making the movement correction directly 
to the signals prior to averaging these. The latter 
method is better in that sense that as the outcome,, a 
signal vector in an illustrative form is obtained. For 
this purpose , minimum norm estimate is used by deter- 
mining the source distribution associated with every 
response to be averaged and by calculating from this a 
signal corresponding to some fixed location of head. 
The problem with the method is the computational slow- 
ness of the minimum norm estimate, the distorting ef- 
fect of the external interferences, and the fact that 
the possibly great distance of the object to the meas- 
urement sensors may distort the result. Methods of 
correcting movements have been described e.g. in the 
publication '^Detecting and Correcting for Head Move- 
ments in Neuromagnetic Measurements", K. Uutela et al, 
Neurolmage, Vol. 14, pp. 1424-1431, 2001. 

The transformation of signal vectors from one 
measurement geometry to another is often necessary 
e.g. when comparing different measurements with one 
another. The question can be e.g. about the transfor- 
mation of each measured signal vector into a signal 
vector of such a measuring device in which the meas- 
urement sensors are situated in determinate positions 
on the surface of some standard object. To calculate 
such virtual signals, the vector field must be divided 
into basis function components, based on which the 
virtual signals can be calculated by forming the basis 
of the vector field for the virtual set of sensors and 
by using estimated components to calculate the virtual 
signals. Conventionally, as the basis, a minimum norm 
estimate basis has been used, in which case the afore- 
mentioned regularisation and interference problems 
hamper the outcome. The use of the minimum norm esti- 
mate for calculating virtual signals has been de- 
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scribed e.g. in the publication ''Trans format ion of 
Multichannel Magnetocardiographic Signals to Standard 
Grid Form", Jussi Numminen et al, IEEE Transactions on 
Biomedical Engineering, Vol. 42, No. 1, 1995. 

The modeling of the sources that produced the 
measured signal is typically implemented by parame- 
trising the source model and by calculating the pa- 
rameterd e.g. by non-linear minimisation of the small- 
est square sum of the error in such a manner that the 
model and the measured signal correspond to each other 
as well as possible. In that case, in the minimisation 
algorithm one must perform the calculation of the pa- 
rameterised model for each measurement sensor specifi- 
cally, which results in a big number of arithmetic op- 
erations when the number of sensors is big. One fur- 
ther problem is the difficulty of setting a good ini- 
tial guesstimate for the model that would fasten the 
convergence of the minimisation algorithm towards the 
correct solution. 

The inaccuracy of the information concerning 
the calibration coefficients and the geometry of the 
measuring device causes errors when analysing an ob- 
ject to be examined. In this sense the measuring de- 
vice must be checked so that the measured signal vec- 
tors are compared to some expected result, in which 
case the calibration coefficients and the geometric 
parameters can be so set that the possible deviation 
from the expected result vanishes. Conventionally, in 
calibration methods one. has used an exactly known sig- 
nal source which is situated in a known position with 
respect to the measuring device. The problem with the 
method is that one must rely on the ideality of the 
signal source. Furthermore, calibration measurements 
of this kind often are cumbersome and time-consuming 
to perform. 
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OBJECTIVE OF THE INVENTION 

The objective of the invention is to elimi- 
nate the drawbacks referred to above, or at least to 
significantly alleviate them. One specific objective 
of the invention - is to disclose a -new type of method 
by which the measured signal can be divided into sig- 
nals associated with the object to be examined and 
into signals that can be. classified as external inter- 
ferences by means of a computationally very easy op- 
eration without the reference sensors. The division of 
the signal .space into these two parts can be performed 
without foreknowledge of the signal sources, so the 
method functions correctly, although e.g. the signal 
sources that are outside the research object and cause 
interferences would change in a non- foreseeable man- 
ner. The signals due to the non- idealities of the 
measuring device can be separated from the interesting 
signals and interference signals because they are 
mainly not included in either one of the aforemen- 
tioned parts of the signal space. In addition, in the 
method in accordance with the invention, there is no 
weakening of the interesting signal nor the illustra- 
tion problems of signals resulting from this, which 
are typical of projection methods. 

Further, the objective of the invention is to 
disclose a faster and more accurate method of compen- 
sation for movement which is used to correct the dis- 
tortion of the signal vector caused by a movement of 
an object to be measured by at the same time correct- 
ing the distortion caused by external interferences 
without any additional calculation. 

Further, the objective of the invention is to 
disclose a calculation method of virtual signals based 
on the fact that in the method, a vector field can be 
presented as components of the spatial series develop- 
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ihent of the field, which components can be limited to 
relate solely to the object to be examined. The compo- 
nents are used to represent the measured signal as a 
signal vector of any sensor geometry. 

Further, the objective of the invention is to 
disclose a method of source modeling in which the 
source distribution of an object to be examined can be 
interpreted by means of the measured basis function 
values, and in some cases it can be solved analyti- 
cally based on these values. 

Further, the objective of the invention is to 
disclose a calibration method of a multi -channel meas- 
uring device that measures an irrotational vector 
field, which calibration method is capable of utilis- 
ing any measured vector field without knowledge of the 
sources that produced the field. In this method, for 
the parameters of the system to be calibrated, such 
values are sought for that the measured signal mainly 
falls in that sub-space of signal space which is 
spanned by the basis functions of the irrotational, 
sourceless vector field. Thus, both the ambiguity as- 
sociated with the sources and the amount of work asso- 
ciated with the preparation of the calibration meas- 
urement are eliminated from the calibration. 

As for the features characteristic of the in- 
vention, reference is made to them in the claims. 

DESCRIPTION OF THE INVENTION 

The invention is based on the utilisation of 
two different sets of vectors in the modeling of a 
signal that originates from two different sources. The 
requirement is that the vector field to be measured is 
irrotational and sourceless. The vectors of the signal 
space in question consist of vector field samples 
taken in different places of the physical space or of 
their linear combinations. In particular, the present 
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invention relates to a method which is used to form 
such a basis in which the information contained in the 
whole vecto^r field can be expressed with a number of 
basis functions as small as possible, in the selection 
of which one has taken into account the irrotational- 
ity and sourcelessness of the vector field. By means 
- of the basis fxinctions in question, the vector field 
can be represented in a form enabling a more enhanced, 
accurate and versatile mathematical operation of the 
field than before in many signal processing and ana- 
lytical applications - 

Furthermore, the invention relates to a 
method the basic idea of .^which is to limit the degrees 
of freedom of the signal vectors of a measuring device 
that measures an irrotational , sourceless vector field 
based on the information on the irrotationality and 
sourcelessness of the field to be measured. In that 
case, the number of possible linearly independent sig- 
nal vectors can be considerably reduced based on the 
limitations set by Maxwell's equations - 

In one embodiment of the present invention, 
the signal to be measured is divided into components 
of two different series developments in such a manner 
that the series developments have different conver- 
gence characteristics in that sense that the signals 
produced by a first source, e.g. an object to be meas- 
ured, and a second source, e.g. external interference 
sources, are included in the convergence areas of dif- 
ferent series developments. In that case, a represen- 
tation is obtained for the signal vectors in which the 
interesting signals and external interferences have 
different basis vectors, the co-ordinates correspond- 
ing to which can be unambiguously solved based on the 
measured signal. 

In one embodiment, the solution is performed 
by minimising e.g. by means of a linear operation the 
difference between a model based on series develop- 
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itients and the measured signal. In case the model de- 
scribes an irrotational and sourceless vector field in 
the area of a measuring device sufficiently perfectly, 
the difference in question is due to a measurement er- 
ror which results solely from the non- idealities of 
the measuring device. Thus, the method can divide the 
signals into three parts: into the interesting signals 
of the irrotational and sourceless vector field, into 
the signals of the irrotational and sourceless vector 
field that can be classified as external interfer- 
ences, and into signals not included in the series de- - 
velopment model and resulting from the non-idealities 
of the measuring device, . such as e.g. noise and 
sources in the measurement area. It must be noted that' 
specifically the claims mention a first and a second 
source, as well as a set of basis vectors, although 
their order is not important from the standpoint of 
the invention. Substantial is that the basis vectors 
describe vector fields originating from different 
sources as accurately as possible. 

Thanks to the taking into account of the 
limitations set by the Maxwell's equations, the number 
of basis vectors needed in the representation of the 
signals is so small that a basis of this kind is prac- 
tical for modern multi-channel measuring devices. In 
addition, the interpretation of the basis vectors is 
physically reasonable and there are a number of solu- 
tions of breaking up signal vectors into basis vectors 
formed in this manner. . The following describes, with 
reference to Fig. 1, the characteristics and utilisa- 
tion of the series developments. It must be noted that 
in all the phases, the limitations set by the Max- 
well's equations are taken into account at the same 
time . 

The present invention has the advantage that 
the method functions as an automatic interference 
elimination method in which the measurement sensors 
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are not separately classified into signal channels and 
reference channels measuring merely interference sig- 
nals, instead all the channels are in the same posi- 
tion when calculating the basis for signal space. The 
automatism of the interference elimination means that 
the method can model external interferences irrespec- 
tive of the structure of the interference field in re- 
lation to the structure of the field produced by the 
object to be examined, or of the variation of the 
structure of the interference field with in respect of 
time, in contrast e.g. to the projection method in 
which the aforementioned signals are projected against 
some predetermined interference sub- space into an or- 
thogonal space. In addition, in the projection method, 
also the signal produced by an object to be examined 
weakens, if it is not orthogonal to the determined in- 
terference sub-space. This is eliminated in the pres- 
ent invention. 

Further, the method in accordance with the 
invention enables one to correct the signal distor- 
tions that are due to the movement of an object to be 
examined faster and more accurately than by previous 
methods, which is a basic requirement of a real-time 
correction for movement. The correction for movement 
is performed by representing the movement of an object 
to be examined as a movement of the measurement sen- 
sors in the vicinity of a stationary object in such a 
manner that the variable geometry is taken into ac- 
count when calculating the basis for signal space. In 
that case, the series development components corre- 
sponding to the object to be examined always are com- 
mensurable irrespective of the position of the object. 
The advantage of the invention is based on . the fact 
that the calculation of the series development compo- 
nents requires quite few computational operations and 
on the fact that the external interferences are mod- 
elled when calculating the values of the series devel- 
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opment components. Thus, this enables one to quickly 
calculate dependable values for the series development 
components corresponding to an interesting source 
area, although the distance and position of the object 
with respect to the measurement sensors would consid- 
erably vary and although the field structure of the 
external interferences would vary during the measure- 
ment . 

Further, as the series development components 
describe the field produced by the sources everywhere 
in the measurement volume, they can be used to calcu- 
late a signal which would be measured by any ficti- 
tious sensor disposed in the same space. Virtual sig- 
nals of this kind are useful e.g. when there is a wish 
to represent the signal produced by an object to be 
examined always in the same geometry irrespective of 
the fact of where the object was really situated dur- 
ing the measurement in relation to the measurement 
sensors. In addition, in virtual signals, the propor- 
tion of measurement errors is smaller than in original 
signals, as in the method one has attenuated the por- 
tion that is not included in the set model of an irro- 
tational, sourceless vector field. Even in the case of 
virtual signals, the invention has the advantage of 
separating external interferences from the signals of 
an object to be examined. 

In addition, the method in accordance with 
the invention enables one to identify individual sen- 
sors the signal measured by which clearly differs from 
the series development model that is based on the sig- 
nals of the whole set of sensors. Thus, the method en- 
ables even a real-time quality control of measurement 
channels and an automatic rejection based on the sig- 
nal analysis. In addition, the co-ordinates of the 
signals with respect to a basis in accordance with the 
invention can be re-calculated after the rejection of 
channels that operate in a non-desired manner, which 



wo 2004/081595 



PCT/FI2004/000144 



13 

enables one to obtain a more accurate result for the 
co-ordinates. 

Further, thanks to the invention, the cali- 
bration of a measuring device or of a measurement ge- 
ometry becomes easier than before. The calibration 
method in accordance with the invention does not need 
a predetermined calibration source, instead it suf- 
fices that the measuring device is known to measure a 
sourceless and irrotational vector field, Advanta- 
geously, the calibration is performed by setting the 
calibration coefficients and geometric parameters of 
the measuring device so that the measured signal vec- 
tors fit in the sub-space spanned by the signal vec- 
tors that are determined by the series developments . 
In other words, the calibration is performed by mini- 
mising the difference between the measured signal vec- 
tors and the sub- space calculated for the signal vec- 
tors that were produced by an irrotational vector 
field, in contrast to the conventional methods which 
minimise the difference between the signal vectors and 
the signal vectors produced by a source that is known 
beforehand as accurately as possible. In the method in 
accordance with the invention, as the measured signals 
it is possible to use e.g. any signals produced by the 
sources in the environment, in which case the calibra- 
tion can be performed completely automatically without 
any preparations of the measurement situation. Fur- 
theannore, the method is more accurate and dependable 
than the previous ones as the inadequacy of the infor- 
mation concerning the geometry of an object to be 
measured does not result in an error of the calibra- 
tion. 

The invention can also be utilised in source 
modeling because the field components calculated irre- 
spective of the external interferences and correspond- 
ing to the target area to be examined represent the 
continuous current distribution of the target area as 
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a finite number of components so that each component 
describes some feature concerning the shape of a cur- 
rent distribution. In . particular, the components can 
be used to check the hypotheses concerning the current 
distribution by comparing the modelled components to 
the theoretically calculated components of the distri- 
butions in - accordance with the hypotheses. In some 
cases, this kind of comparison is easier than when us- 
ing as the measurement result, the measured signal 
vector values, as it is done in several conventional 
methods. In some cases, it is possible to analytically 
calculate from the series development components the 
geometry of the current distribution without non- 
linear minimisation algorithms, which makes the deter- 
mination of the distribution clearly faster than in 
the previous methods. 

The invention also provides a new way of de- 
signing the geometry or set of sensors of a measuring 
device by optimising the capability thereof to de- 
pendably measure the series development components of 
an interesting object up to an order as high as possi- 
ble by maximising at the same time the difference of 
the signal spaces associated with the object of the 
measuring device in question and the external inter- 
ferences, such as e.g. the angle between the spaces. 
In this manner it is possible to obtain a set of sen- 
sors which is capable of measuring even the most com- 
plicated signals in respect of their field shapes as 
accurately as possible, while at the same time sepa- 
rating the interesting signals from the external in- 
terferences. 

LIST OF FIGURES 

In the following, the invention will be de- 
scribed by means of detailed examples with reference 
to the accompanying drawing, in which 
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Fig. 1 shows one generalised embodiment of 
the present invention; and 

Fig. 2 is a schematic representation illus- 
trating one embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

At first we describe the mathematical back- 
ground of the invention with reference to Fig. 1. To 
illustrate the idea presented in this invention, in 
Fig. 1, the space has been schematically divided into 
five areas by means of concentric spherical shells. 
The sensors (si, s2, s3, s4, s5, s6) of the multi- 
channel measuring device are situated in area 3, which 
is left in between the second and third innermost 
spheres. The sensors are connected to an analysis de- 
vice 10. In area 2 there is a current loop I^^ which 
represents an interesting source of magnetic field 
situated in the measurement object. A current loop 
situated in area I^^ represents a source that produces 
external magnetic interferences. 

In the following, the representations of the 
field produced by these two sources are examined by 
such series developments the development point of 
which is the common midpoint of the spheres inside 
area 1. It is commonly known that the field produced 
by source I^^ can be represented as a Taylor's develop- 
ment which is convergent in area 1, which thus is such 
an area each point of which is closer to the develop- 
ment point than the closest source point. The field 
produced by source .1^^ can also be represented as a 
convergent, asymptotic series development outside area 
2, i.e. in an area in which all the points are farther 
from the development . point than any source point. 

An asymptotic series development is herein 
used to mean e.g. a series development that is based 
on spherical harmonics, the terms of which are proper- 
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tional to the negative integer - powers r" of the dis- 
tance r calculated from the development point. The 
same remarks also apply to a magnetic field produced 
by an external interference source I^^ : . the Taylor's 
series development thereof is convergent everywhere in 
areas 1, 2 and 3, whereas the asymptomatic series de- 
velopment is convergent in area 5. 

As can be seen from Fig. i, in area 3 con- 
taining the sensors (si, . s2, s3, s4, s5, s6) , where 
the field is thus measured, the convergent series de- 
velopment of an interesting signal is asymptotic, 
whereas the convergent series development of an exter- 
nal interference in this area is the Taylor's series. 
When for the measured signal (signal space), both a 
Taylor's basis and a basis based on an asymptomatic 
series development are created, the signals coming 
from an interesting object can be separated from in- 
terferences coming from external sources based on the 
fact that the external interferences are mainly shown 
in the sub- space spanned by the basis vectors that are 
based on the terms of the Taylor's series, whereas the 
interesting signals are included in the sub-space de- 
termined by the basis vectors associated with the 
asymptomatic series. 

In the following, one advantageous 
embodiment of the invention is shown with reference to 
Fig. 2. Let us think about a situation in which from 
the space, a point is chosen in such a manner that 
the smallest and biggest distance of the measurement 
points from this are r„,„ ja r^, respectively. 

The measurement points are designated in Fig. 
2 as points on a curved surface. Let us assume that 
there are in the space two sources of an irrotational 
vector field, the first of which is situated at a dis- 
tance r, < r^^ from point r„ and the latter at a dis- 
tance > from point ^. 
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There, is some series development T which is 
convergent inside such a sphere whose radius is the 
distance from the development point to the" closest 
singular point, as well some other series development 
M which is convergent in an area which encompasses a 
space containing singular points. Thus, the vector 
field produced by source 1 can be expressed in the 
area of the measurement points as a development of the 
M series with respect to point r^, but not as a devel- 
opment of the T series, and correspondingly, the vec- 
tor field produced by source 2 can be expressed in the 
area of the measurement points as a T development with 
respect, to point but not as an M development. Thus, 

the signal vectors produced by sources 1 and 2 can be 
separated from one another, if the measured signal 
vector is presented in a basis whose basis functions 
correspond to the independent components of the T and 
M development with respect to the chosen development 
point. It must be noted that the development point 
does not need to be the same for different series de- 
velopments and that there are several different alter- 
natives for the developments. It suffices that the 
convergence areas of the developments are the same as 
with the aforementioned T and M developments. The 
method described becomes general in a situation in 
which there are several sources both at distance r^^ < 
and at distance r^^ > r^. The development point can 
be placed e.g. so that all the sources associated with 
the object to be examined are situated in area r^^ < 
and the external interferences in area r , > r . Fig. 2 
shows two spheres which describe the convergence areas 
of the series developments. Inside the smaller sphere, 
the T developments of the vector fields produced both 
by source 1 and source 2 are convergent ♦ with respect 
to point ( designated with Tl, T2) . In the space be- 
tween the smaller and the bigger sphere, the M devel- 
opment of the vector field produced by source 1 is 
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convergent, and the T development of the vector field 
produced by source 2 is convergent (designated with 
Ml, T2) . Outside the bigger sphere, the M developments 
of the vector fields produced by both of the sources 
are convergent with respect to point ( designated 

with Ml, M2) . The fields produced by the sources can 

be separated based on the series developments in the 

space between the spheres, where also the measurement 

points are situated. Thus, when the signals of the 

measurement points are expressed at the same time both 

as the T and M series developments with respect to 
point To , source 1 solely affects the signal through 

the components of the M development, and source 2 
solely through the components of the T development. 

If the aforementioned series development com- 
ponents form a linearly independent basis for signal 
space, then the interesting signals can be separated 
from the interferences independently. To reduce the 
effect of noise, one may have to use regularisation, 
in case the susceptibility to interference of the ba- 
sis is big. Herein the susceptibility to interference 
is used to mean the product of the norms of the matrix 
and of its reverse matrix. The linear independence and 
the susceptibility to interference of the basis depend 
on the geometry of the measurement arrangement and on 
the number of the components to be modelled in rela- 
tion to the number of measurement points. The inven- 
tion provides a basis for the signal vectors of an ir- 
rotational vector field that is more practical than 
before and in which the regularisation means substan- 
tially the selection of the number of series develop- 
ment components to be included. This is easy to solve 
because of the fast attenuation of the components as 
the order of the development increases. 

The following describes the division of the 
measured signal vectors into an object to be examined 
and into external interferences using as the example 
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the Taylor's series development as well as the multi- 
pole development of the spherical harmonics. As the 
vector field being measured functions the magnetic 
flux density. Let it be assumed that the development 
point, for the series developments . can be selected in 
such a manner that the sources associated with the ob- 
ject-to be measured are situated inside an r ' radial 
sphere, and the sources associated with the external 
interferences are situated outside an R radial sphere 
when the origin of both of the spheres is the develop- 
ment point in question. If r < and R > r^, when 
and are correspondingly the smallest and the big- 

gest distance of the development point to the measure- 
ment channels, then the division into interesting sig- 
nals and external interference signals can be per- 
formed as shown in the following. 

The Taylor •s series development is convergent 
inside an R radial sphere when R is the distance from 
the development point of the series to the closest 
singular point, herein to a source. Thus, the signal 
vector ip^^j produced by a source s^^ situated at a 

distance r*^^ © R can be expressed with respect to the 
development point by means of the Taylor series of the 
magnetic field, but the Taylor development is not con- 
vergent for a signal vector produced by a source s^^ 
situated at a distance • r. Instead, for ^„ it is 

possible to use e.g. the multi-pole development. The 
signal measured on channel k is the surface integral 
of a normal component of the magnetic flux density 
over the measurement area. Now we consider the fre- 
quency without limiting the planar sensors, in which 
case the surface integral in question can be written 
in the following form: 



wo 2004/081595 



PCT/FI2004/000144 



20 



(l>,=n^jb{r)dS, (1) 

Sk 

The co-ordinates of the" Cartesian co- 
ordinates are designated with symbols x^, x^, x,. In 
that case, the Taylor development of the component 
of the magnetic flux density at the measurement point 



b (n = f f f i£245^^ 

:^&opi mlllpldx'^dyldz' ' 



(2) , 



wherein the partial derivative concerns the 
development point of the series development. Thus, 



wherein 



G = 



^dx By Bz 



dbydbydby 

ydx dy dz 



^ db.db^db, 
^ dx dy dz 



nz 



nz 



(nx+ ra/+ nz) 



(3) 



(4) 



and 



(5) 



Thus 

S/c 



(6) 



The integral expression of the previous equa- 
tion is a vector, so it can be denoted as 
^outk = ^c* Gpfc 

(7), 

in which case the signal vector produced by 
an external source is 
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lout=dLog(N^GP)^ (Qj^ 
wherein 

N =^[n^..jtf^] and P = [Pi -Pw]- 

The. Maxwell's equations Vxfc=0 and, V*^ = 0 
that are valid in the measurement area limit the num- 
ber of independent G matrixes in such a maimer that in 
the Taylor development of nth order there are 
n7.=(n + 2)^-1 of independent terms, in which case each 

signal vector produced by external sources can be ex- 
pressed as follows : 

.Tip 

lout = ^ Xj^diogiN^GiP) = ^ A^Wt 

1=1 i=l , (9) 

The basis formed by the signal vectors S„ is 

denoted by matrix 3^, in which case the signal vectors 
produced by external signal sources are of the form: 

^ouf=SrXr. (10) 



In a mult i -pole development, the magnetic 
scalar potential is expressed by means of so-called 
even and odd spherical harmonics. On the other hand, 
the magnetic fliox density is the negative gradient of 



the potential in question, so herein 
pressed in the form: 



can be ex- 



1=0 m=0 



1 



H 



(11) 



wherein a,, and. b,, are mult i -pole coefficients, 
and the even and odd spherical harmonics can be ex- 
pressed by means of Legendre's polynome P^^icoa*) in 
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the form (6, 9) = P,^co^)co^mp) 

Yim (e,<p) = Pjj;co^)sir(mp) 

The angles " and 'P are angles of bearing of 
the spherical co-ordinates. The redundancy due to the 
irrotationality and sourcelessness of the field is 
automatically noticed in a multi-pole development, and 
the dimension of the nth order multi-pole basis is 
nw=(n + 1)^-1 when the mono-pole term (n = 0) is not no- 
ticed. Now the basis vectors ^ corresponding to the 
multi-pole development can be formed by the same prin- 
cipe as in the case of the Taylor basis, whereby a 
multi-pole basis- is obtained. Thus, the signal vec- 
tors associated with an object to be examined are of 
the form : 

^^"^^"^ (12) 



The actual main point of the invention is the 
combining of the bases ja S„ into one basis for sig- 
nal space in such a manner that 5 = [.S,„5r] and jc = 



whereby the measured signal vector can be represented 
in the form: 

This equation describes a signal vector meas- 
ured by an ideal measuring device from an irrotational 
and sourceless vector field. In practice, a non- ideal 
measuring device produces into the model of an irrota- 
tional and sourceless vector field an inaudible meas- 
urement error ®, so by practical measuring devices the 
foregoing equation has to be written in the form: 



(14) 
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The series development components can be 
solved by minimising the difference between the meas- 
ured signal and the model by the method of the small- 
est sc[uare sum, which in this case corresponds to the 
multiplying of the measured signal vector by the 
-pseudo inverse of the matrix S: 



The solution of an ideal case is unambiguous, 

in case the dimension of the measured signal vector ^ 
is bigger than the number of vector fields of the ba- 
sis S and in case the basis S is linearly independent. 
These requirements are typically realised by present- 
day multi-channel measuring devices, if the limita- 
tions set by the Maxwell's equations are taken into 
account. Furthermore, to ensure an error-free esti- 
mate, a sufficiently big number of series development 
components must be included in the basis S. As these 
requirements are fulfilled, the measured signal vector 
can always be divided into the part produced by an ob- 
ject to be measured, into the part produced by exter- 
nal interferences, and into the part not included in 
the irrotational vector field (measurement error) : 



Ls^Sr^, (17) 
c = 0-Sjc (18) 



The field of application of the invention is 
wide as in the method, from the measured signal vec- 
tor, a natural representation is calculated containing 
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substantially the same information as the signal ved- 
tor, but in a form independent of the sensor geometry 
and typically in a lower- dimensional form than the 
measured signal vectors. Furthermore, the components 
associated with an object to be examined and with . the 
external interferences are separate, unlike in the 
signal vector, which is the sum of an interesting sig- 
nal and of interferences. 

In the following, one implementation mode of . 
the invention is described by means of an example. Let 
it be assumed that we wish to measure the magnetic 
field produced by the brain of a testee who does not 
keep his or head still, instead constantly moves in 
the vicinity of the measuring device. Furtheirmore , in- 
terference sources external of the testae's head are 
associated with the measurement situation, the struc- 
ture of the magnetic field produced by which changes 
as a function of time. Interference sources of this 
kind can be e.g. the magnetic particles in a moving 
testee ' s clothes and the testee • s heart . 

The magnetic field measurements of the brain 
are usually made as repetition measurements in which 
the response produced by the brain and relating to the 
desired stimulus is averaged several times in a row to 
improve the signal-to-noise ratio. In the case of the 
testee described above, the averaged signal vector to 
be obtained as the outcome is the mean value of the 
signal vectors produced by a source that was situated 
in different places and is physiologically invariable 
from one stimulus to another, and thus the average 
signal cannot be used in the analysis as such. Fur- 
thermore, the elimination of external interferences by 
a projection method based on a predetermined interfer- 
ence stib-space is out of the question, since the field 
structure of the interferences has varied in a non- 
foreseeable manner. An interference elimination based 
on separate reference sensors is not justified either. 
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since the testee can time to time be far from both the 
actual signal and reference sensors, whereby a divi- 
sion- into an interesting signal and interferences that 
is based on shear uniformity of the field or on a sig- 
nal amplitude is not very successful. Instead, the in- 
terference elimination can be performed by the series 
development method .described above by averaging the 
component estimates associated with an interesting ob- 
ject, in the calculation of which, one has taken into 
account the place and position of the testee with re- 
spect to the measuring device at different moments of 
time. 

. .The motion of a testee is modelled in such a 
manner that the motion of the head is represented as a 
motion of a set of sensors in the environment of an 
immovable head. In that case, at a moment of time t 
basis matrixes S^(t) and S^(t) are calculated, from 
which a basis S(t) is formed, based on which it is 
possible to calculate, by the pseudo inverse, an esti- 
mate x{t) in accordance with equation 15. When the co- 
ordinates of the measuring device are connected to the 
co-ordinates of a stationary head, the portion Jc,„(r) of 

the component estimate x(t) associated with the object 
is independent of the testee »s motion. Thus, instead 
of the averaging of the signal vectors associated with 
the stimuli it is worth averaging the estimates 

whose mean value contains the same information as the 
mean value of the signal vectors, but without the dis- 
tortion caused by the testee 's motion. The outcome of 
a repetition measurement is thus a performed correc- 
tion for movement and an automatic elimination of ex- 
ternal interferences, so the mean value can be used as 
such for the analysis without other pre-processing op- 
erations . 

One analysis method includes the use of 
multi-pole components in source modeling, which is 
particularly useful when the source distribution asso- 
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ciated with an object to be examined can be solved 
analytically from the estimated components. This is 
successful at least in a case in which the field is 
produced by one current dipole inside a spherical vol- 
ume conductor. On the other hand, the outcome can also 
be represented as a virtual signal vector of any sen- 
sor geometiry by calculating the basis of the sensor 
geometry in question, whereby the virtual signal vec- 
tor is obtained in accordance with equation 16 by us- 
ing an averaged estimate jc„ . 

As in the calculation of the component esti- 
mates one must use information on the geometry and 
calibration coefficients of the measuring device, it 
is important that the geometric parameters and the 
calibration coefficients are known as accurately as 
possible. In the method of the invention, these are 
found out by solely utilising the information on the 
fact that the measuring device measures an irrota- 
tional vector field. In that case, each measured sig- 
nal vector shall be included in the sub- space S, and 
systematic deviations from this condition refer to an 
incomplete piece of information on the calibration and 
geometry of the measuring device. Thus, the calibra- 
tion can be performed by setting the calibration coef- 
ficients and the geometric parameters so that the sub- 
space S calculated based on these corresponds to the 
sub-space spanned by the measured signal vectors. 

The invention is not limited merely to the 
examples of its embodiments, instead many variations 
are possible within the scope of the inventive idea 
defined by the claims. 
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CLAIMS 

1. A method for analysing a magnetic signal 
and/or field registered by a multi -channel measuring 
device in a geometry in which an interesting source is 
in measurement volume vi, the sensors measuring the 
field or the components thereof outside volume VI in- 
side volume 'V2,-and the sources of magnetic interfer- 
ences outside volume VI and V2 in volume V3, which can 
be infinite, characterised in that 

- parametrising in volume V2 the magnetic 
field produced by the interesting sources disposed in 
volume VI as a sum of elementary fields, each of which 
is. irrotational,. sourceless and infinite outside vol- 
ume VI in such a manner that for the parametrised mag- 
netic field, a representation of the desired accuracy 
in volume V2 is obtained; and 

- parametrising in volume V2 the magnetic 
field produced by the interference sources disposed in 
volume V3 as a sum of elementary fields, each of which 
is irrotational, sourceless and infinite outside vol- 
ume VI in such a manner that for the parametrised 
field, a representation of the desired accuracy in 
volume V2 is obtained; and 

- calculating the signal vector of the meas- 
uring device corresponding to each elementary field; 
and 

- measuring the magnetic signal by sensors; 

and 

separating the fields that were produced 
from sources disposed in different volumes by calcu- 
lating the components of the measured signal vector in 
the basis formed by the signal vectors associated with 
the elementary fields. 

2. The method as defined in claim 1, 
characterised in that the signal vector is 
formed as a sum of the desired components of the meas- 
ured signal vector. 
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3. The method as defined in claim 1, 
characterised in that the signal vector is 
formed as a sum of the desired components of the meas- 
ured signal vector in a set of sensors which is in a 
position deviating from its position during the meas- 
urement with respect to the interesting source. 

4. The method as defined in claim 1, 
characterised in that the signal vector is 
formed as a sum of the desired components of the meas- 
ured signal vector in such a set of sensors which de- 
viates from the set of sensors used in the measurement 
with respect of its geometry. 

5. The method as defined in claim 1, 
characterised in that as the elementary 
fields associated with the sources disposed in volume 
VI and volume V3 , the fields associated with the terms 
of the series development that represents the harmonic 
potential of two fields that are convergent in the 
measurement space on different conditions are used, 
whereby the magnetic signals originating from volume 
VI and volume V3 are separated as sets of basis vec- 
tors of their own. 

6. The method as defined in claim 5, 
characterised in that the series develop- 
ment estimating the field potential generated from the 
sources disposed in each volume VI and V3 is selected 
so that it is convergent in volume V2 as fast as pos- 
sible, in which case a number of basis vectors as 
small as possible suffices to obtain a description as 
accurate as possible for the signals caused by differ- 
ent, sources , 

7. The method as defined in claim 5, 
cha.racterised in that as the series devel- 
opment estimating the field potential generated from 
the sources disposed in each volume VI and V3, either 
the Taylor's series, the series based on spherical 
harmonics or some other applicable series is chosen. 
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8. The method as defined in claim 1 or 1, 
characterised in that in the selection of 
the elementary fields, the rediindancy due to the har- 
monics of the potential is taken into account. 

9. The method as defined in claim 5, 
characterised in that the development point 
and/or . points of the series developments are chosen 
from volume VI or V3 . 

10. The method as defined in . claim 1, 
characterised in that the signal corre- 
sponding to the signal source disposed in volume VI or 
the signal corresponding to the signal source disposed 
jLn volume VS. is separated from the .measurement signal 
as a function of time. 

11. The method as defined in claim 1, 
characterised in that the components of the 
measured signal vector in the basis formed by the sig- 
nal vectors associated with the elementary fields are 
deducted from the measurement signal, in which case 
the separation signal corresponds to the signal origi- 
nating from volume V2 . 

12. The method as defined in claim 11, 
characterised in that the signal originat- 
ing from volume V2 is registered as a function of 
time. 

13. The method as defined in claim 11, 
characterised in that a channel that func- 
tions in an unusual manner is identified by means of 
the separation signal. 

14. The method as defined in claim 13, 
characterised in that the identified chan- 
nel functioning in an unusual manner is taken into ac- 
count in the further analysis in accordance with the 
invention with a weighting coefficient different with 
respect to the other channels. 

15. The method as defined in claim 1, 
characterised in that the motion of the 
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source is perceived, and based on this information, 
the elementary fields associated with the source are 
represented in the co-ordinates of the source, 

16 • The method as defined in claim 1, 
characterised in that the interesting 
sources are in volume V3 and the interference sources 
in volume VI . 

17. A method for calibrating a multi-channel 
measuring device that measures an irrotational, sour- 
celess vector field by means of measurement signals, 
characterised in that 

- selecting a set of basis vectors corre- 
sponding to irrotational, sourceless elementary fields 
in measurement volume V2 ; and 

setting the calibration coefficients and 
geometric parameters of the measuring device to values 
which minimise the deviation of the signal vector from 
the signal space spanned by a basis in accordance with 
the set of basis vectors. 

18. The method as defined in claim 15, 
characterised in that the change is per- 
ceived by monitoring the place of a known signal 
source arranged in a first and/or second space with 
respect to the measurement geometry. 

19. A device for analysing a magnetic signal 
and/or field registered by a multi-chaxmel measurement 
device, the device comprising: 

a set of sensors (si, s2, . . .sn) ; and 
an analysis device (10) to which the sensors 
are connected, in which case the measurement is ar- 
ranged in a geometry in which the interesting source 
is in measurement volume VI; the sensors measuring the 
field, or the components thereof outside volume VI in 
volume V2; and the sources of the magnetic interfer- 
ences outside volume VI and V2 in volume V3, which can 
be infinite, characterised in that the 
analysis device (10) is arranged to 
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- parametrise in volume V2 the magnetic field produced 
by the interesting sources disposed in volume VI as a 
sum of elementary fields, each of- which is- irrota- 
tional, sourceless ajid finite outside volume VI in 
such a manner that for the parametrised magnetic 
field, a representation of the desired accuracy in 
— volume V2 is obtained; and 

- parametrise in volume V2 the magnetic field 
produced by the interference sources disposed in vol- 
ume V3 as a sum of elementary fields, each of which is 
irrotational, sourceless and finite outside volume V3 
in such a manner that for the parametrised field, a 
representation of the desired accuracy in volume V2 is 
obtained; 

- calculate the signal vector of the measur- 
ing device corresponding to each elementary field; and 

- measure the magnetic signal by sensors; and 

- separate the fields that were produced from 
sources disposed in different volumes by calculating 
the components of the measured signal vector in the 
basis formed by the signal vectors associated with the 
elementary fields. 

20. The device as defined in claim 19, 
characterised in that the signal vector is 
formed as a sum of the desired components of the meas- 
ured signal vector, 

21. The device as defined in claim 19, 
characterised in that the signal vector is 
formed as a sum of the desired components of the meas- 
ured signal vector in a set of sensors which is in a 
position deviating from its position during the meas- 
urement with respect to the interesting source. 

22. The device as defined in claim 19, 
characterised in that the signal vector is 
formed as a sum of the desired components of the meas- 
ured signal vector in such a set of sensors which de- 
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viates from the se.t of sensors used in the measurement 
with respect to its geometry. 

23. The device . as defined in claim 19., 
characterised in that as the elementary 
fields associated with the sources disposed in volume 
VI and volume V3 , the fields associated with the terms 
of the series development that represents the harmonic 
potential of two fields that are convergent in the 
measurement space on different conditions are used. 

24. The device as defined in claim 23, 
characterised in that the series develop- 
ment estimating the field potential generated from the 
sources disposed in each volume VI and V3 is selected 
so that it is convergent in volume V2 as fast as pos- 
sible, in which case a number of basis vectors as 
small as possible suffices to obtain a description as 
accurate as possible for the signals caused by differ- 
ent sources. 

25. The device as defined in claim 23, 
characterised in that as the series devel- 
opment estimating the field potential generated from 
the sources disposed in each volume VI and V3, either 
the Taylor's series, the series based on spherical 
harmonics or some other applicable series is chosen. 

26. The device as defined in claim 19 or 25, 
characterised in that in the selection of 
the elementary fields, the redundancy due to the har- 
monics of the potential is taken into account. 

27. The device as defined in claim 23, 
characterised in that the development point 
and/or points of the series developments are chosen 
from volume VI or V3 . 

28. The device as defined in claim 19, 
characterised in that the signal corre- 
sponding to the signal source disposed in volume VI or 
the signal corresponding to the signal source disposed 
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in volume V3 is separated from the measurement signal 
as a function of time. 

29. The device as defined in claim 19, 
characterised in that the components of the 
measured signal vector in the basis formed by the sig- 
nal vectors associated with the elementary fields are 
deducted from the measurement signal, in which case- 
the separation signal corresponds to the signal origi- 
nating from volume V2 . 

30. The device as defined in claim 29, 
characterised in that the signal originat- 
ing from volume V2 is registered as a function of 
t ime - . 

31. The device as defined in claim 29, 
characterised in that a channel that func- 
tions in an \inusual manner is identified by means of 
the separation signal, 

32. The device as defined in claim 31, 
characterised in that the identified chan- 
nel functioning in an unusual manner is taken into ac- 
count in the further analysis in accordance with the 
invention with a weighting coefficient different with 
respect to the other channels. 

33. The device as defined in claim 1, 
characterised in that the motion of the 
source is perceived, and based on this information, 
the elementary fields associated with the source are 
represented in the co-ordinates of the source. 

34. The device as defined in claim 19, 
characterised in that the interesting 
sources are in volume V3 and the interference sources 
in volume VI . 

35. The device as defined in claim 19, 
characterised in that 

selecting a set of basis vectors corre- 
sponding to irrotational , sourceless elementairy fields 
in measurement volume V2; and 
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- setting the calibration coefficients and 
geometric parameters of the measuring device to values 
which minimise the deviation of the signal vector from 
the signal space spanned by a basis in accordance with 
the set of basis vectors - 

36 • The method as defined in claim 33, 
characterised in that the change is per- 
ceived by monitoring the place of a known signal 
source arranged in a first and/or second space with 
respect to the measurement geometry. 

37. The device as defined in claim 19, 
characterised in that just the desired com- 
ponents of the measured signal vector are saved. 

38. The device as defined in claim 19, 
characterised in that in locating an inter- 
esting source, just the desired components of the 
measured signals are used. 
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